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We present results for the effect of a second massive quark in the relation between the 
pole and the minimal subtracted quark mass at the three loop level. We also consider 
the analogous effect for the wave function renormalisation constant. Some technical 
details of the calculation are given. Our result is phenomenologically relevant for the 
bottom quark including virtual charm effects. 

1 Introduction 

Quark masses are fundamental parameters of the Standard Model (SM) and thus it is de- 
sirable to determine their numerical values with the highest possible precision. In order to 
do so it is necessary to fix a renormalisation scheme which defines the quark mass. Often 
physical observables are expressed through some "short distance" mass [I] to obtain predic- 
tions which are free of the renormalon problem. To do so. one frequently needs the relation 
between the on-shell and the MS mass. Many authors contributed to the latter task [2l[3l|4]. 
In this contribution wc present the recently published calculation [5] , where we have included 
the effect of a second nonzero quark mass to the relation between the quarks in those two 
schemes at the three loop level. After having reduced the problem to the calculation of mas- 
ter integrals we use two different ways to solve them, namely the differential equation and 
the Mellin-Barnes method. From the phenomenological point of view this result is important 
for the bottom-quark including effects from virtual charm-quarks. As a byproduct we also 
obtain the corresponding contribution to the wave function renormalisation constant. 

2 Renormalisation constants 

Introducing the decomposition of the quark self-energy 

T,{q,mq) = mgl]i(q^,m,) + - mq)Y,2{q^ ,mq) , (1) 
we can express the renormalisation constants, which are defined through 



by [im 



TO,,o = Mq , = V^V- , (2) 



= l + ^,{MlMq), (3) 
^2°')"' = l + ^M^qTr^^iiq^Mq) ^ +j:2{MlMq). (4) 



dq 



Tp is the quark field renormalised in the on-shell scheme with mass m^, Mq is the on-shell 
mass and bare quantities are denoted by a subscript 0. Thus, to obtain Z^^ one only needs 
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Figure 1: Master integrals. The solid/wavy lines correspond to particles with mass Mq/0, 
the dashed lines denote the quark with the second nonzero mass. 



to calculate Si for ~ M^. To calculate Z2^, one has to compute the first derivative 
of the self-energy diagrams. The mass renormalisation is taken into account iteratively by 
calculating one- and two-loop diagrams with zero-momentum insertions. 

In the case of the mass renormalisation it is convenient to consider the ratio between the 
on-shell and MS renormalisation constants 

which is finite. Here we have introduced the renormalisation dependent MS-mass mq(^). 

In contrast to the wave function renormalisation constant contains next to ultra- 
violet also infrared divergences. Thus it is not possible to construct a finite quantity by 
considering the ratio between the on-shcU and MS renormalisation constant. 

3 Computational techniques 

In order to compute the on shell self energy we use QGRAF \7\ to generate the feynman 
diagrams and the various topologies are identified with the help of q2e and exp [SI [9]. In a 
next step the reduction of the various functions to so-called master integrals (Mi's) has to 
be achieved. For this step we use the so-called Laporta method [lOl [TT] which reduces the 
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three-loop integrals to 26 MPs. We use the implementation of Laporta's algorithm in the 
program Crusher [12| . It is written in C++ and uses GiNaC |13] for simple manipulations 
like taking derivatives of polynomial quantities. In the practical implementation of the 
Laporta algorithm one of the most time-consuming operations is the simplification of the 
coefficients appearing in front of the individual integrals. This task is performed with the 
help of Fermat [14] where a special interface has been used (see Ref. |T5]). The main 
features of the implementation are the automated generation of the integration-by-parts 
(IBP) identities [16] . a complete symmetrisation of the diagrams and the possibility to use 
multiprocessor environments. 

In Fig. [T] a graphical representation of the master integrals can be found. As indicated 
in the figure, many topologies contain more than one master integral. We have chosen two 
independent ways to compute the e-expansion of the master integrals. The first one relies 
on the differential equation method [17]. With this method we were able to evaluate all 
but four master integrals in analytic form. With the help of our second method, based on 
the Mellin-Barnes technique (see, e.g., Ref. [18]) we can get numerical results for all master 
integrals. Here we have used the Mathematica package MB.m |19j . 



3.1 DifFerential-Equation-Method 

First we set up differential equations in z = Alf/Mq, where Mf is the second nonzero 
quark mass, for each of the individual integrals. Each equation will contain the master 
integral itself and integrals belonging to the corresponding topology. The latter can again 
be reduced to the set of master integrals. For all the topologies which contain only one 
master integral (e.g. all six-propagator integrals in Fig. [ij this gives a "simple" equation 
whereas for the other topologies wc get two or three coupled equations, respectively. In 
the next step we expand the differential equations in e. Choosing an appropriate basis for 
the master integrals, all equations decouple order by order in e. We can now solve for all 
integrals by repeated use of Euler's variation of the constant method. The initial conditions 
we need are all known from the literature [20l [211 IHl [6] . 

We were able to get analytical results for all master integrals in terms of Harmonic 
Polylogarithms (HPL's [23]) up to order e^^ . Unfortunately we were not yet able to get 
analytical results in higher orders in e for the four integrals belonging to the two topologies 
depicted as the leftmost ones in the second row of Fig. [TJ For all other integrals we provide 
analytical results up to the order we need in our calculation. To calculate the MPs in terms 
of HPL's it is necessary that the (pseudo)poles in the corresponding differential equations 
are all of the form 1/z, l/(z -I- 1) and l/(z — 1), were these poles can occur up to arbitrary 
order. This is the case in all topologies except for the ones mentioned above. We have not 
found a proper change of variables to bring the differential equation into this form. As a 
consequence wc only managed to integrate the integrals in these topologies up to the order 
e~^. The remaining parts can e.g. be integrated numerically with Mathematica. 

To evaluate our results numerically and for general algebraic manipulations of terms 
involving HPL's we use the Mathematica package HPL [24l [25]. 
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3.2 Mellin-Barnes-Method 

The McUin-Barncs method as a tool for the evaluation of Feynman integrals has become 
very popular in the recent years. The basic formula is |18| 



which transforms a propagator like term into a complex contour integral. A common recipe 
to evaluate Feynman integrals is the following: First one introduces Feynman parameters for 
a loop variable. Then one can perform the corresponding momentum integration. After that 
one applies formula ([6]) to the denominators containing the Feynman parameters. Finally 
the Feynman parameters can be integrated yielding the Mellin-Barnes representation of the 
original integral. This procedure has recently been automatised [26]. 

The Mellin-Barnes integration is to be performed along a contour which reaches from 
— oo to oo on the imaginary axis with indentations such that the poles of r(. . . + z) and 
those of r(. . . — z) are separated by the contour. 

Mellin-Barnes integrals usually have singularities for certain values of their parameters. 
If there are for example Gamma functions of the form r(£ -I- z) T{—z), it is not possible to 
find an appropriate integration contour when e ^ 0. The integral is therefore singular in e 
and this is how UV poles of Feynman integrals manifest themselves in their Mellin-Barnes 
representation. One has thus to regularise the integral, that is, identify the e poles. This 
can be done by shifting the integration contour using the residue theorem. Prescriptions 
to do so have been given in Refs. [27] and [28], the latter has been automatised, see Refs. 



Finally the regularised integrals can be expanded in e and evaluated by numerical inte- 
gration, which is also implemented in the package MB .m |19| . or by application of the residue 
theorem and summing up the residues. Depending on the dimension of the integrals and 
the complexity of the integrand this can be done numerically or analytically. 

To calculate the master integrals for this work we first simplified the Mellin-Barnes 
integrals by inserting the representations of known subtopologies. The regularisation has 
been done partly using MB . m and the approach of Ref . [27] . One- and two-dimensional MB- 
integrals were calculated via their sum representation, higher dimensional integrals using 
MB.m. 

The 4-line integrals can all be represented as Mellin-Barnes-integrals of maximal dimen- 
sion 1, which can be evaluated as single sums. For the 5 line master integrals we find 
representations of dimension 2 to 5. The integrals with 6 lines have 3 to 5 dimensional 
representations. We find good agreement with the results obtained from the differential 
equation method. 

Inserting the results for the master integrals into the final result we observe large numer- 
ical cancellations near Mf = between the contributions originating from different master 
integrals. On the other hand, the expansion for Mf/Mq <^ 1 converges very fast, which is 
relevant for Mf = rric and Mq = mt,. For this reason we decided to derive an expansion of 
our result including eighth order terms. The coefhcients that could not be obtained analyt- 
ically can be evaluated numerically from their Mellin-Barnes-representation, which is in our 
case at most two-dimensional. 





[291 dl]. 
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4 Results and applications 



As an application of our result wc want to discuss the charm quark effects in the relations 
between the pole, the MS and the IS quark mass. For illustration we use mb{mb) = 4.2 GeV, 
nT-cirric} = 1-3 GeV, /i = mi, and ai^''(mb) = 0.2247. The relation between the on-shcU and 
the MS mass leads to 



4.2 + 0.401 



0.199 + 0.0094 



0.145 + 0.0182 



GeV, (7) 



where the tree- level, one-, two- and three- loop results are shown separately. The contribu- 
tions from the charm quark mass which vanish for rric are marked by a subscript nic. 
One observes that the higher order contributions are significant. In particular, the two-loop 
charm quark effects amount to 9 MeV and the three-loop ones to 18 McV. A similar bad 
convergence is observed in the relation between the IS" mass [30] and the pole mass 

Mb. For Mb = 4.8 GeV, mc(mc) = 1.3 GeV, = Mb and ai^\Mb) = 0.2150 it is given by 



M, 



IS 



4.8 - 0.049- 0.073 + 0.0041 



- 0.098 + 0.0112 



GeV. 



(8) 



However, the relation between the 15* and the MS quark mass is much better behaved as can 
be seen in the following example where we have chosen M^^ = 4.69 GeV, m,c{mc) = 1.3 GeV, 
(Mi-^) = 0.2167 



M^^ and ai*'^' 



nib = 



4.69 - 0.382 - 0.098 + 0.0047 



- 0.030 + 0.0051 



GeV. 



(9) 



The two-loop charm effects amount to only 4.7 MeV and three-loop ones to 5.1 MeV. We 
want to mention that in case only the linear approximation [3j of the charm quark mass 
effects is used the corresponding three-loop results in Eqs. ^ and ([9]) read 0.0167 and 
0.0037, respectively. 
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